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Difference control schemes for controlling unstable fixed points become important if the exact 
position of the fixed point is unavailable or moving due to drifting parameters. We propose a 
memory difference control method for stabilization of a priori unknown unstable fixed points by 
introducing a memory term. If the amplitude of the control applied in the previous time step is 
added to the present control signal, fixed points with arbitrary Ljapunov numbers can be controlled. 
This method is also extended to compensate arbitrary time steps of measurement delay. We show 
that our method stabilizes orbits of the Chua circuit where ordinary difference control fails. 

PACS numbers: 05.45.+b, 84.30.Ng, 07.50.Ek 



INTRODUCTION 



II. STABILIZATION OF FIXED POINTS BY 
DIFFERENCE CONTROL 



Two main classes of control strategies are established 
in chaos control: The OGY control algorithm almost 
a standard method for controlling chaos, does not pro- 
vide any readjustment of the fixed point during in-situ 
measurements without loss of control. In many experi- 
mental systems, however, it is desirable to use a control 
strategy that does not rely on the knowledge of the ex- 
act position of the fixed point, because the location of 
the fixed points can change due to drifts in parameters. 
On the other hand, the time-continuous control method 
proposed by Pyragas is practically limited by the re- 
quired sampling rate, and does not allow stabilization of 
arbitrary orbits as has recently been shown in §j . 

Both problems are circumvented by simple time- 
discrete difference control 0. It is limited to a certain 
range of Ljapunov numbers. Control of arbitrary peri- 
odic orbits can be achieved if the algorithm is applied 
only every second period 0,(3 or by a memory method. 



In this paper we present an improved memory differ- 
ence control (MDC) method that takes control ampli- 
tudes into account that were applied at previous time 
steps. MDC allows one to stabilize drifting fixed points 
with arbitrary Ljapunov numbers and shows an enlarged 
region of stability. 

This method is generalized when dealing with measure- 
ments delayed by r time steps (orbit periods). This task 
is accomplished by increasing the number of memorized 
control amplitudes by r. Given the stable and unsta- 
ble directions of the fixed point with sufficient accuracy, 
only one accessible control parameter for each unstable 
direction is required to achieve control. 

We compare difference control and MDC at the well- 
known Chua oscillator and show that orbits for which 
difference control fails are stabilized by MDC. 



In experimental situations a Poincare section is com- 
monly used to reduce the dynamics to a time-discrete 
description by an iterated map 



xt+i = f{xt,r). 



(1) 



Here r denotes a set of control parameters that are in 
the unperturbed dynamics assumed to be constant or 
varying on a slow time scale compared to the the length 
of a period orbit. 

The idea to control chaos by small perturbations 
of control parameters implies that r becomes time- 
dependent. The time-dependent control parameter rt 
is updated at each discrete time step defined by the 
Poincare section. Its value is determined according to 
the specific control algorithm and is kept constant for a 
part of the orbit. Without loss of generality we choose 
X* — Q and r = for the fixed point to be stabilized. The 
linearized equation of motion around the unstable fixed 
point then becomes 



Xt+i = Lxt -f Mrt 



where 
L^j := 



dxj 



and Mi- 



S*=0,f=0 



(2) 



(3) 



The Ljapunov numbers of the orbit are the eigenvalues 
of L. Here one has to distinguish the Ljapunov num- 
ber of an orbit (or time-discrete map) from the local 
(or conditional) Ljapunov exponent and the commonly 
used global Ljapunov exponent being an ergodic aver- 
age over the attractor In principle it is possible to 
proceed with a multiparameter control by using as many 
control parameters as there are degrees of freedom, i. e. 
dim(7^) = dim(a;). Instead, it is common to follow Ott, 



2 



Grebogi and Yorke to transform the system to the 
eigensystem of L. Control is applied only in the unstable 
subspace SJ . The evolution of the equation of motion is 
again of the form of eq. Q with reduced dimension of 
L. 

In difference control 0| the control parameter is up- 
dated at the Poincare section according to 



n = K{xt - xt^i). 



(4) 



In contrast to OGY control, difference control is limited 
to fixed points with Ljapunov numbers between —3 and 
— 1 0, ■ A stability diagram 5, ^ for the case of one 
unstable eigenvalue is shown in Fig. ^ 




FIG. 1: Stability area for difference control (one unstable 
Ljapunov number L): For \L\ < 1 the system is stable without 
control. Fixed points with Ljapunov number —3 < L < —1 
can be controlled it M ■ K is chosen within the area bounded 
by the triangle. The line within the triangle shows the optimal 
value for M ■ K. 

Simple difference control uses information that is par- 
tially out of date, resulting in an additional degree of 
freedom from the delayed amplitude Xt+i- This fact is 
illustrated by imagining two agents controlling the same 
system in turns. If they do not communicate, control 
fails because of the inherent delay of the system to be 
controlled. This effect is compensated when using the 
information contained in the previous value of the con- 
trol amplitude rt_i. 



III. MEMORY DIFFERENCE CONTROL 

We define memory difference control by |in| | 

ft = K{xt - xt-i) -f Nrt-i. 



(5) 



Combined with l^l, we obtain a dynamical system which 
reads in delayed coordinates for x and r 



n 




-MK MN 


-K N 




In order to stabilize the fixed point all eigenvalues ai 
of the matrix in Q must have a modulus smaller than 
one. This ensures exponential convergence to the fixed 
point. If all parameter values are chosen such that all 
eigenvalues become zero, the fixed point is reached after 



a finite number of time steps. In fact this can be guaran- 
teed by MDC. We first assume that M and (L — 1) are 
both invertible, and that the number of accessible control 
parameters is equal to or greater than the dimension of 
the unstable manifold, i. e. dim(r) > dim(a5'). Then all 
eigenvalues are zero [ll|, if 



K ^ -M-^L^{L- 1)^1 
TV = M-^L{L - 1)-^M. 



(7) 



The concept of MDC can be generalized to stabiliza- 
tion of (known and unknown) fixed points when x can be 
measured only after a finite number of delay steps [Tllll3j| : 
If the system is measured with r steps delay, (Q is re- 
placed by 



Vi 



<i<^ 



K 

Nr+l 



~-M-^L^+^{L-l)-^ (8) 
-Nr^VM (9) 
M-^L^+^{L~l)-^M (10) 



where the feedback now contains a sum over r 
preceding control amplitudes: 



1 



T+l 



= K{xt-r - Xt-r~l) + ^ NiXt^r^ 



(11) 



A similar control scheme can be applied for OGY control 
by choosing K = -M-^L^+i and Ni = ~M-^L'M, 
(1 < I < t). For details see [TllIT^ . 



IV. THE STABILITY DIAGRAM IN THE 
ONE-DIMENSIONAL CASE 

Since the optimal control values are never exactly 
matched in experiments, it is important to know the 
complete stability diagram in {K,N), in particular the 
optimal value of K for given N and vice versa. In the 
one-dimensional case, the characteristic equation is given 
by 

= a[a^ -{L + MK + N)a + [MK + NL)]. (12) 

The stability region in the (if, A^)-plane, i.e. where all 
eigenvalues have modulus smaller than one, is the trian- 
gle shown in Fig. [3 Its corners are given by 



{MK, A^)+i,+i 
{MK, A^)-i,-i 
(Mif,iV)+i,_i 



(1-L,1) 

{L + lf (L + 3) 



(L-l)'(L-l) 
(-1-L,1) 



(13) 
(14) 
(15) 



where the eigenvalues take the values -1-1 and —1 as in- 
dicated by the indices. Two sides of the triangle are de- 
termined by the conditions that one eigenvalue is equal 



3 



to +1 and —1, respectively. The third is given by 
MK + NL = +1 where the eigenvalues are a complex 
conjugate pair on the unit circle. The line that deter- 
mines that N with minimal eigenvalues for a given K 
(and vice versa) is given by the algebraic expression 

{MK f + 2MKN + TV^ + {2L - A)MK - 2LN + 



VCR- 



Ci 



where both eigenvalues coincide. 
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FIG. 2: Stability region of memory difference control for 
L = -2 in the {MK, iV)-plane. Within the triangle MDC 
is stable; the contour lines show the decrease of the largest 
eigenvalue to zero in {MK, N) = (4/3, 2/3). 



V. THE CHUA OSCILLATOR 

We demonstrate the efficiency of the improved dif- 
ference method by controlling unstable periodic orbits 
of the chaotic attractor of Chua's oscillator 0. The 
Chua circuit is an autonomous system. The Poincare 
section necessary for control is obtained by an electronic 
zero-crossing-detector. The frequency is in the range of 
V = 2t: l^/LCi ~ 3.6kHz and allows the usage of digital 
signal processing tools to implement control algorithms. 

Parameter drifts, e.g. temperature drifts, naturally oc- 
cur in electronic circuits and difference control methods 
have the advantage that they follow the drifting fixed 
point. In order to get access to an appropriate control 
parameter, a VCR (voltage controlled resistor) has been 
added to the circuit. The basic dynamics is nevertheless 
determined by the resistor R. 

Furthermore the Chua circuit allows to investigate in- 
teresting ranges of Ljapunov numbers simply by choosing 
different values of the main control parameter R. 

The dynamics of the Chua circuit is described, in first 



FIG. 3: The Chua circuit: The negative resistance (NIC) 
is both nonlinearity and energy source of the circuit. Rough 
adjustment of the control parameter can be done by adjust- 
ing R. Control is applied with the voltage-dependent resistor 
(VCR). 

approximation, by the normalized equations 

ii — a{v — M — f{u)) 

V — u — V + w (16) 

w = —I3v 

where / is the input-output function of the negative re- 
sistance approximately described by the piecewise linear 
descending function 

f{u) = mow + ^(mi - mo)(|u + 1| - |u - 1|) (17) 

with Too > TOi [l^ . Rather than solving these equations 
numerically, we demonstrate the stabilization of an un- 
stable periodic orbit directly in the electronic system. 




0.5 1 1.5 2 

x(t) 

FIG. 4: Return map Xt+i{xt) in the Poincare surface of sec- 
tion for three different values of the control parameter. From 
this data one obtains values of L and M to adjust the optimal 
parameters of control. 



VI. IMPROVED DIFFERENCE CONTROL OF 
THE CHUA SYSTEM 

The standard control strategy is to measure the re- 
quired system variables, generate the Poincare map for 
e. g. three adjacent values of the control parameter 
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(Fig. Q , and to calculate the parameters of the feedback 
to the control parameter (Fig.[Sl). 

In the present case the return map xt+i{xt^ zt) is ap- 
proximately a function of xt alone. Therefore only two 
variables are required. The first one, j/t, is used to deter- 
mine the Poincare surface of section by a zero crossing 
detector. The second one, xt, is used for the calculation 
of the control. 



chaotic circuit 



Poincare section 



control 



X 



DSP system 



FIG. 5: Implementation of the control method: One coordi- 
nate is used to generate the Poincare section. A digital signal 
processor is used for measurement, computation, and appli- 
cation of the control amplitude. 
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FIG. 7: Stabilization of an orbit in the double-scroll chaotic 
attractor: In agreement with theory, control without memory 
(N=0) fails due to a Ljapunov number L < —3. Within the 
triangle the orbit is stabilized by memory difference control. 
Again, attractor and orbit are shown in the inset. 



The stability region for different values of the mem- 
ory term TV and feedback gain K is measured by chang- 
ing the values until control is entirely lost. The lower 
bound of K is easily recognized by a doubling of the sta- 
bilized period. However, the upper bound of K, where 
the loss of control is noise-induced, is more difficult to es- 
timate. In Fig.lHlthc stability region for a stabilized orbit 
in the single-scroll chaos is shown. The corresponding 
Ljapunov number L = —2.15 ± 0.04 has been calculated 
from the Poincare map. The stability diagram includes 
the stability region of simple difference control as the 
special case {N = 0). Stabilization of the same periodic 
orbit in the double-scroll chaotic attractor is not possible 
with simple difference control. This is due to a Ljapunov 
number of L = —3.27 ± 0.08 for which the method is 
predicted to fail. 




4 6 
feedback gain K 

FIG. 6: Stabilization of an orbit in the single-scroll chaotic 
attractor: Within the measured triangle memory difference 
control stabilizes the orbit of the Chua circuit. The special 
case of simple difference control is given by A'' = 0. The inset 
shows the attractor and the stabilized periodic orbit. 

The stability region of the stabilized orbit in the single- 
scroll (Fig. O and double-scroll attractor (Fig. O have a 



broad overlap. Thus it is possible to choose universal 
values of {K, N) that allow tracking of an orbit from one 
regime to the other without changing parameters of the 
controlling circuit. 




3 4 5 6 
feedback gain K 

FIG. 8: Measured Ljapunov exponents of the control tran- 
sients for N — and N = 0.7 (single scroll regime) compared 
to theoretical Ljapunov exponents (lines) from the eigenvalue 
analysis. 



The significant improvement of memory difference con- 
trol compared to simple difference control is demon- 
strated by the estimation of Ljapunov exponents (con- 
traction rates) from the transient. Fig. 8 shows the sta- 
bility regimes for different feedback gains K. Simple 
difference control corresponds to iV = 0, and MDC to 
N — 0.7. The range of controllability is broadened and 
the Ljapunov exponents are smaller, equivalent to faster 
convergence. The measurements are in good agreement 
with our theoretical predictions (Fig. |SJ). Due to noise- 
induced loss of control, it was impossible to obtain reli- 
able measurements for large K of the case N = 0.7. 



5 



VII. RE-ESTIMATION OF LJAPUNOV 
NUMBERS FROM THE BORDERS OF 
STABILITY 

From the borders and corners of the measured stabil- 
ity region the exact values of the Ljapunov number of 
the controlled cycle are re-estimated, similar to the ap- 
proach used in [l^ for OGY control of the Henon map. 
Since two values of N are exactly one, four of the six co- 
ordinate values can be used to determine L and M by a 
least-square fit weighted by the variances of the measured 
values giving 

Lss = -2.069 ± 0.03 Lds = -3.24 ± 0.03 

Mss = 0.376 ± 0.015 A/ds = 0.488 ± 0.005 

(18) 

for the orbit stabilized in the single scroll (ss) and double 
scroll (ds) attractor, respectively. These values are in 
good agreement with the values given in Section lVll which 
were obtained from the Poincare map (Fig.^. 

VIII. CONCLUSIONS 

The dynamical behavior and stability conditions of dif- 
ference control and memory difference control are fun- 
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